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Abstract 

■ Using the normally ordered Gaussian form of displaced-squeezed thermal field characteristic 

of average photon number n, we introduce the photon-added squeezed thermo state (PASTS) 
and investigate its statistical properties, such as Mandel's Q-parameter, number distribution (as 
a Legendre polynomial), the Wigner function. We then study its decoherence in a photon- loss 
channel in term of the negativity of WF by deriving the analytical expression of WF for PASTS. 
It is found that the WF with single photon-added is always partial negative for the arbitrary 
' values of n and the squeezing parameter r. 
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^ ! 1 Introduction 



Nonclassicality of fields has been a topic of great interest in quantum optics and quantum information 
processing [1 . Experimentally, the traditional quantum states, such as Fock states and coherent 
states as well as squeezed states, have been generated but there are some limitations in using them 
for various tasks of quantum information process [2]. Alternately, it is possible to generate and 
manipulate various nonclassical optical fields by quantum superpositions and subtracting or adding 
photons from/to traditional quantum states 0|ll3IS[71llllS[Ii[ni[Tl[T3l[Tl[T51[TB]. 

On the other hand, the single mode displaced, squeezed, mixed Gaussian states have been paid 
enough attention by both experimentalists and theoreticians. Marian et. al [T7l[T8j investegated the 
superposition of a squeezed thermal radiation and a coherent one. They examined the squeezing 
properties of the field using the distribution functions of the quadratures. For Gaussian squeezed 
states of light, a scheme is also presented experimentally to measure its squeezing, purity and 
entanglement [TH1 HO]- As is well known, dissipative quantum channels tend to deteriorate the 
degree of nonclassicality (i.e., render quantum features unobservable) . Thus, it is usually necessary 
to investigate the decoherence properties in dissipative channels, such as dynamical behaviors of the 
partial negativity of Wigner function (WF) and how long a nonclassical field preserves its partial 
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negativity of WF. For instances, the nonclassicality of single photon-added thermal states in the 
thermal channel is investigated by exploring the volume of the negative part of the WF [21] ; Souza 
and Nemes have derived an upper limit for the mixedness of the single bosonic mode Gaussian states 
[2"2] in a thermal channel. 

In this paper, we shall introduce the photon-added squeezed thermo state (PASTS) and investigate 
its statistical properties, such as Mandel's Q-parameter, number distribution, the Wigner function. 
We then study its decoherence in a photon-loss channel in term of the negativity of WF by deriving 
the analytical expression of WF for PASTS. It is found that the WF with single photon-added is 
always partial negative for the arbitrary values of n and the squeezing parameter r. The work is 
arranged as follows: In section 2, we introduce the state PASTS and derive its normalized constant, 
and photon number distribution is discussed in section 3. Section 4 is devoted to calculating the 
WF. In the last section, we explore the decoherence of PASTS in a photon-loss channel by discussing 
the evolution of WF. 



2 PASTS and its normalization 

For a displaced-squeezed thermal field, the density operator is 

Ps = D((3)S{r)p c S\r)D\pl (1) 

where D(f3) = exp(/3a^ — (3*a), and S(r) = exp[ir(QP + PQ)/2], are the displacement operator and 
the squeezing operator respectively, (3 = (q + ip)/V2, Q = s ^-, P = [a, a f ] = 1, and 



p c = (l-e <^)e 

(k is the Boltzmann constant, T denoting temperature), is qualified to be a density operator of 
thermal (chaotic) field, since tvp c = 1. For a coherent state \z) = exp [z(V — z*a) |0) [251 due to 
a\z) — z\z), matrix elements of any normally ordered operators : O (a' , aj : (the symbol : : denotes 
normally ordering) in the coherent state is easily obtained, i.e, 



{z\:6{a\a) : \z') = O (z* , z') {z \z') 

= 0{z*,z')exp\ ' +z*z'\, (2) 




so in Ref. [27[ 128] by using the Weyl ordering invariance under similarity transformations and the 
technique of integration within an ordered product of operators (IWOP) Fan et al have converted 
p s to its normally ordered Gaussian form 



P> = — ' CX P \ ^2 t 1 > ( 3 ) 



l_ \ {q-Qf {p-Pf 

nr 2 - eXP \ 2r x 2 2r| 

where 

2tI = (2n + l)e 2r + 1, 2r| = (2n + l)e- 2r + 1, (4) 

and n is the average photon number for p c , i.e. n = ( e huJ / kT — i)- 1 ^5]. The form in Eq.Q is 
similar to the bivariate normal distribution in statistics, which is useful for us to further derive the 
marginal distributions of p s . 

Theoretically, the PASTS can be obtained by repeatedly operating the photon creation operator 
on a displacement squeezed thermal state, so its density operator is defined as 
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Pm = N- l a! m Ps a m , (5) 

where m is a non-negative integer, N m =tr(a' m p s a m ) is the normalization constant. Using Eq.((3]) 
we known immediately the normally ordered Gaussian form of p m , i.e., 

" - N - a^i-^-^l)a-:. (6) 



tit 2 r \ 2t{ 2t. 
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Next we shall determine the normalization constant N m . Using the completness relation of 
coherent states J \z) (z\ = 1 as well as Eq.©, we have 



trpn 



, f d 2 2 

tr( Pm J — 



d 2 z 



TlT 2 J TT 



N- 



z\ : a) m exp 



(q-Q) 2 (p-Pf 



2rl 



2rf 



TlT 2 J TT 



-z* m z m exp 



where we have set 



-A \z\ z + B*z + Bz* + Cz 2 + Cz* 2 + D 



1 ( <1 ip 



A 2T 2 + lr B - V2{tI + rl>- 



C 



Due to trp m = 1, thus we know 



1 1 

4^f + 4^f' 



D 



2 2 



N m = 



T\T 2 



TlT 2 8 A 

Further using the following integral formula 



J ^_L z *rn z va ^ |^|2 + B * z + Bz * + Q -,2 + (J z *2 j 

J ^ exp j-A \z\ 2 + B*z + Bz* + Cz 2 + Cz* 2 } 



{-l) m e D d m 



f d 2 z ( ~i 1 

/ exp C \z\ +£,z + V z* + fz 2 + gz* 2 = exp 

J n I J v/C 2 - Afg 



C 2 - 



(7) 



(8) 



(9) 



(10) 



whose convergent condition is Re(C ± / ± g) < and Re ( ^±f±g ^ < 0, Eq.|9]) can be rewritten as 
follows 



N„ 



(A 2 - 4C 2 y 1/2 exp 



A\B\* + CB* 2 + CB 
A 2 - AC 2 



(11) 



which is the normalization constant of PASTS for photon-added number m. In particular, when 
(3 = leading to B = D = 0, Eq. tfTTj) reduces to the following form, 



TiT 2 8A m V ' 



Especially when m = 0, 1, 2 then the normalization constants are given by Nq = 1, iVi 
and N 2 = I (3 Tl 4 + 2rf r| + 3t|) , respectively. 



(12) 



i W + T 2 2 ) 
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To see clearly the photon statistical properties of the PASTS, we will examine the Mandel's 
Q-parameter defined as 

(a) 2 a 2 ) 

Qm = Vtt^ - ( flta ) > ( 13 ) 

(aJa) 

which measures the deviation of the variance of the photon number distribution of the field state 
under consideration from the Poissonian distribution of the coherent state. If Qm = we say the 
field has Poissonian photon statistics while for Q > 0(Q < 0) we say that the field has super- 
(sub-) Poissonian photon statistics. From Eq. (fTT|) and N m ~ti(a' m p s a m ), we can easily calculate 
( ta) = ^±1 - 1, and (a^a 2 ) = - ^±1 + 2;t hus we obtain the Q-parameter of the PASTS 

„ N m+2 -AN m+1 +2N m N m+1 -N m 

Qm = N N N ' ( ' 

It is well known that the negativity of the QM-parameter refers to sub-Possonian statistics of 
the state. But a state can be nonclassical even though Qm is positive. This case is true for the 
present state. From Fig.l, one can clearly see that for the cases of m = (Fig. 1(a)), Qm is always 
positive; while for m/0 (for instance m — 1) and a given n value, Qm becomes positive only when 
the squeezing parameter r is more than a certain threshold value that increases as m increases. In 
addition, from Fig. 1(a) and Fig. 1(b) one can see that the threshold value of r decreases as n increases. 
We emphasize that the WF has negative region for all r and n thus the PASTS is nonclassical (see 
next section below). 

3 Photon number distribution of PASTS 

In this section, we study photon number distribution of PASTS optical field. Using the un- normalized 
coherent state |a) = exp[aat] |0), leading to \n) = J^^ n l a ) l a=0 > (( a l a ') = e ° " ) ; ^ is easy to 
see that the photon number distribution formula is given by 



V(n) = U(p m \ri) (n|) = (n\p m \ri) 
1 1 d^ n 

= 1 j * nA ,» {a\p m \a!) \ a=a '=o ■ (15) 

t\T2 n\ da "da " 

Employing the normal ordering form of p m in Eq.@, Eq. (|15p can be put into the following form 

KT-lpD j2n 

= J T T da *n da ,n W* m a' m exp [B*a' + Bat + Ca' 2 + Cot 2 + (1 - A) a* a'] } | a=Q , =0 . 

(16) 

Further expanding the exponential term exp [(1 — A) a* a'] as series and using the generating function 
of single-variable Hermite polynomials, 

Hn(x) = exp (2xt - t 2 ) | t=0 , (17) 
we can calculate the photon number distribution (PND) of PASTS, i.e., 
_ N^e y> (1 - A) 1 d m+l d m+l 

1=0 

,2h 



exp{B*a' + Ba* + Ca 12 + Ca* 2 } \ a=a , =0 



da'™ da 

N-i\C\ n e D ^{l-Ay 



1 z 1=0 
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dB m+l H n 



±B/(2VC) 



2 



(18) 
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After making the scale transformation and noticing the recurrence relation -^;H n (x) = ^™}y H n ^i(x), 
we can easily obtain 



V(n) 



N- l e D n!(l - A) 1 |C| r 



T\T 2 



E 

J=0 



Z![(n-m-Z)!] 



m—l 



LB/(2VC) 



Especially when /3 = 0, Eq. lflT)]) reduces to 



(19) 



V{n) 



(_-]\m V [(«-m)/2] n!(l-A)'"-'"- 2 ^|C|^ 

^(AZ-4C*r 1/2 
( n _ m )!^_(yl2_4 C .2)-i/2' 



where 



o = y/(l-A) 2 -4C* = ^/(r?-!)^ 2 - l)/(nr 2 ), 



and in the last step of ([20]) we have used the new expression of Legendre polynomials [8] 

[™/2] l\l 



E 

Z=0 



m i(i _ 

2 2l {l\) 2 (m-2l)V 



(20) 



(21) 



(22) 



Eq. ([Hi)) or (|2H)) is just the the analytical expression of the PND of PASTS. In particular, when 
m = 0, Eq.® becomes (with = 0) 



rir 2 a 



(23) 



Eq. ([23|) is just the PND of the squeezed thermo state which seems a new result. The PNDs of 
PASTS for some given parameters (n, r) and m are plotted in Fig. 2. From Fig. 2 it is found that 
the PND is constrained by n m. By adding photons, we have been able to move the peak from 
zero photons to nonzero photons (see Fig. 2 (a)-(c)). The position of peak depends on how many 
photons are created and how much the state is squeezed initially. In addition, comparing Fig. 2(b) 
and Fig. 2(d) we see that, for a given m, the "tail" of PND becomes more "wide" with the increasing 
parameter r. 



4 Wigner function of PASTS 

The Wigner function (WF) [3T] was first introduced by Wigner in 1932 to calculate quantum cor- 
rection to a classic distribution function of a quantum-mechanical system. It now becomes a very 
popular tool to study the nonclassical properties of quantum states. It is well known that WFs 
are quasiprobability distributions because it may be negative in phase space }32j . Nevertheless, the 
partial negativity of the WF is indeed a good indication of the highly nonclassical character of the 
state. Thus, to study the dynamical behaviors of the partial negativity of WF and understand that 
a nonclassical field preserves its partial negativity, Wigner distribution may be very desirable for 
experimentally quantifying the variation of nonclassicality [33 . 

The presence of negativity of the WF for an optical field is a signature of its nonclassicality. In 
this section, using the normally ordered form of PASTS, we evaluate its WF. For a single-mode 
system, the WF in the coherent state representation \z) is given by pH] 



■5 



W(a,a*) 



a 2\ a \ 2 



-z\pm \z)e 



-2(za* —z*a) 



(24) 



where a = (x + iy) /y2. Then substituting Eq. ((6]) into Eq. (|24l) and using Eq. ([2]), we derive the 
WF of PASTS 



W(a,a*) 



TXT\T2 OF" 



(^2_4C 2 )" 1/2 exp 



-F\E\ 2 + E* 2 C + E 2 C 
F 2 - AC 2 



where we have set 



F = 2 - A, E = B-2a 
and used Eq. fTO]) . Especially when (3 = 0, Eg. ([25]) reduces to 



(-l) m e 2 H aF^ 



{< 



W(a,a*) 

VI 

and further when to = 0, 1 Eq. (|27[) becomes to 

1 



exp 



-4F|g| 2 +4ct C+4a 2 C 
F a -4C a 



]} 



7^7 (^ 2 - 4C2) 



-1/2 



W,„ =0 (a,a 



7r(2n+ 1) 

which is the WF of the squeezed thermo state, and 



exp 



e 2r x 2 + e 2r y 2 
2n + 1 



W m=1 (a, a ) = exp 



e 2r x 2 + e 2r y 2 
2n + 1 



(25) 
(26) 

(27) 

(28) 
(29) 



respectively, where we have set 



1 2r 1 2 e- 2 '' , . 

M= T ^ 30 

(2rl+l) 3 t 2 + t 2 

91 = 31 

(2n+l) 3 t 2 + t 2 

1 T 2 -I- T 2 - 4t 2 T 2 

T = 1 / 1+ : 2 I 12 - (32) 

(2n+l) 3 t 2 +t 2 

Eq. (]25|) just agrees with the result of Eq.(48) in Ref. [27], whose form is normal distribution. 

From Eq. (f27]) one can see that the WF of the PASTS is always real, as expected. When the factor 
Sola: 2 + my 2 + T < in Eq.®, the WF of the PASTS with to = 1 has its negative distribution 
in phase space. Noticing 971, 91 always positive, this indicates that the WF of the PASTS always 
has the negative values under the condition T < (i.e., (t 2 + t| — 4r 2 r| < 0)) at the phase space 
center q = p = 0. In fact, by substituting Eqs. ((4]), (]30]) - (132p into T < 0, we find that for the 
arbitrary values of n and r, the WF with m = 1 is always partial negative. 

Using Eq.f27]). the WFs of the PASTS are dipicted in phase space for several different values of 
to, n and r in phase space. Fig. 3 exhibits the WFs of the PASTS in phase space with to = 1 for 
different n, r. It is easy to see that the WFs of single-PASTS always have the negative region. The 
minimum in the negative region becomes larger with the increasing of n (see Fig3.(a) and (c)). In 
Fig. 4, we have presented the WFs with n = 0.5, r = 0.3 for different m, which indicates that the 
peak (absolute) value of WF become smaller as the increasing parameter m. The partial negativity 
of WF indicates the nonclassical nature of the PASTS field. 
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5 Evolution of WF in a photon-loss channel 



When the PASTS evolves in the amplitude decay channel, the evolution of the density matrix can 
be described by the following master equation in the interaction picture |_35j , 



dp 



K (2a pa) — a) ap — pa) a) , 



(33) 



where k represents the rate of decay. By using the thermal field dynamics theory and thermal 
entangled state representation, the time evolution of WF at time t to be given by the following form 
[36] . i.e., 



W(a,t) = - —e-^ 



T 

T=l-e 



-2Kt 



W(z,0). 



(34) 



Eq. (f3"4")) is just the evolution formula of WF of single mode quantum state in photon-loss channel. By 
observing Eq. (pS4")) . we see that whent — > 0,T — > 0, ^ cxp(-|. \a - ze~ Kt \ 2 ) — >• 8 (a - z) 5 (a* - z*) , 
so W (a,t) — » W (a,0) as expected. Thus the WF at any time can be obtained by performing the 
integration when the initial WF is known. 

For simplicity, here we only discuss the special case (3 — 0. Substituting Eq. ([27|) into Eq. ([34|) and 
using Eq. (flQ|) , we derive the time evolution of WF for PASTS in photon- loss channel: 



W(a,t) = 



2(-l) m OF* 



ttT 



(35) 



where we have set 



N 



F 2 - AC 2 



4F + hs^l _ 2 ) (F 2 - AC 2 ) 



64C 2 



4Ne- 



r 2 

16NCe- 2Kt 
T 2 



AF 



2e 



T 



[F - AC 2 



2 



(36) 



In particular, when m = Eq. (|35[) becomes 

(A 2 - AC 2 ) 



W m=0 (a, t) = 



J|a| 2 +ka* 2 +ka 2 



ttT 



(37) 



which is just the WF of the squeezed thermo state in photon-loss channel. This result can also be 
checked by substituting Eq.flgS]) into Eq. ([3"4]) . 

When nt exceeds a threshold value, the WF has no chance to be negative in the whole phase 
space. At long time nt — > oo, leading to N —> ^ p _ 2 ^_ 16C , 2 > K — > —2, k — > 0, the WF in Eq. 
becomes 



e -2|«| 2 

W(a,oo) = , (38) 

which corresponds to the Gaussian state. In Fig. 5, the WFs of PASTS are depicted in phase space 
with to = 1 and r = 0.3 for several different nt. It is easily seen that the negative region of WF 
gradually disappears as nt increases. This implies that the system state reduces to a Gaussian state 
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after a long time interaction in the channel. Thus the loss of channel causes the absence of the 
partial negativity of the WF if the decay time kI exceeds a threshold value. 

In Figs. 6, we have also presented the time-evolution of WF for different r. One can see clearly 
that the partial negativity of WF decreases gradually as r increases. The squeezing effect in one of 
the quadratures can be seen in Fig. 6. In Eq. (|35|) , for the PASTS we have obtained the expression 
of the time evolution of WF. In principle, by differentiating as shown in Eq. (|35[) we can derive 
the WF of other PASTS (m 1). But its form is too complicated. However, we can draw the 
Wigner distributions of PASTS for m = 2,3, 5 by numerical simulation, as shown in Fig. 7 (a)-(c), 
respectively, from which one can see that the absolute value of the negative minimum of the WF 
decreases as m increases. 

6 Conclusions 

Based on the normally ordered Gaussian form of displaced-squeezed thermal field, we have introduced 
a kind state: the photon-added squeezed thermo state (PASTS). Then we have investigated the 
statistical properties of PASTS (such as Mandel's Q-parameter, number distribution, the Wigner 
function) and its decoherence in photon-loss channel with dissipative coefficient n in term of the 
negativity of WF by deriving the analytical expression of WF for PASTS. It is found that the 
photon number distribution is just a Legendre polynomial and that the WF with single photon- 
added is always partial negative for the arbitrary values of n and the squeezing parameter r. The 
technique of integration within an ordered product of operators brings convenience in our derivation. 
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Figures caption: 

Fig.l Mandel's Q-parameter of PASTS as a fuction of r with m — 0,1,5,10,30 (from top to 
bottom) for (a) n = 0.3; (b) n = 1. 

Fig. 2 Photon number distributions of PASTS with n = 1 for (a) r = 0.3, m = 0; (b) r = 0.3, m = 
l;(c) r = 0.3, m = 5;(d) r = 0.8, m = 1. 

Fig. 3 WF of PASTS for m = 1 (a) n = 0, r = 0.3; (b) n = 0,r = 0.8;(c) n = 0.5, r = 0.3;(d) 
n = 0.5, r = 0.8. 

Fig. 4 WF of PASTS for n = 0.5, r = 0.3 (a)m = 0; (b) m = 2;(c)m = 3 ;(d) m = 5. 
Fig. 5 The time evolution of WF of PASTS for m = l,r = 0.3, and n = 0.5 with (a) Kb = 0.05; 
(b) Kt = 0.15;(c) nt = 0.2;(d) Kb = 0.4. 

Fig. 6 WF of PASTS for m = 1, n = 0.5 at nt = 0.05, with (a) r = 0.01; (b) r = 0.5;(c) r = 1. 
Fig. 7 WF of PASTS for r = 0.3, n = 0.5 at Kb = 0.05, with (a)m = 2; (b) m = 3. (c) m = 5. 
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Figure 2: Photon number distributions of PASTS with n = 1 for (a) r = 0.3, m — 0; (b) r — 0.3, m : 
l;(c) r = 0.3, m = 5;(d) r = 0.8, m = 1. 
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Figure 3: WF of PASTS for m = 1 (a) ft = 0, r = 0.3; (b) ft = 0,r = 0.8;(c) ft = 0.5, r = 0.3;(d) 
n = 0.5,r = 0.8. 




Figure 4: WF of PASTS for ft = 0.5, r = 0.3 (a)m = 0; (b) m = 2;(c)m = 3 ;(d) m = 5. 




Figure 5: The time evolution of WF of PASTS for m = 1, r = 0.3, and n = 0.5 with (a) /rf = 0.05; 
(b) Kt = 0.15;(c) Kt = 0.2; (d) Kb = 0.4. 
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Figure 6: WF of PASTS for m = 1, n = 0.5 at nt = 0.05, with (a) r = 0.01; (b) r = 0.5;(c) r = 1. 
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Figure 7: WF of PASTS for r = 0.3, n = 0.5 at nt = 0.05, with (a)m = 2; (b) m = 3. (c) m = 5. 
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